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ABSTRACT 

Infrared dark clouds (IRDCs) are cold and dense reservoirs of gas potentially available to form stars. Many of these clouds are likely to be pristine 
structures representing the initial conditions for star formation. The study presented here aims to construct and analyze accurate column density 
and dust temperature maps of IRDCs by using the first Herschel data from the Hi-GAL galactic plane survey. These fundamental quantities, are 
essential for understanding processes such as fragmentation in the early stages of the formation of stars in molecular clouds. We have developed 
a simple pixel-by-pixel SED fitting method, which accounts for the background emission. By fitting a grey-body function at each position, we 
recover the spatial variations in both the dust column density and temperature within the IRDCs. This method is applied to a sample of 22 IRDCs 
exhibiting a range of angular sizes and peak column densities. Our analysis shows that the dust temperature decreases significantly within IRDCs, 
from background temperatures of 20-30 K to minimum temperatures of 8-15 K within the clouds, showing that dense molecular clouds are not 
isothermal. Temperature gradients have most likely an important impact on the fragmentation of IRDCs. Local temperature minima are strongly 
correlated with column density peaks, which in a few cases reach N// 2 = 1 x 10 23 cm -2 , identifying these clouds as candidate massive prestellar 
cores. Applying this technique to the full Hi-GAL data set will provide important constraints on the fragmentation and thermal properties of 
IRDCs, and help identify hundreds of massive prestellar core candidates. 
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1. Introduction 

One of the major issues in star formation is understanding the 
physical conditions, including the temperature and density struc- 
ture, of the material out of which stars will form. While fun- 
damental, these quantities are difficult to determine. These re- 
gions are cold and therefore radiate energy away through dust 
at the far infrared/submillimeter wavelengths, a frequency range 
barely accessible from the ground because of the opacity of the 
atmosphere. 

For the first time, a high angular resolution view of the 
Universe from 70 to 500 //m is now becoming available thanks 
to the ESA Herschel Space Observatory ( [Pilbratt & et al.|2010| ). 
The fre quency range provided by its two photometry instruments 
PACS ( |Poglitsch&etal.|2010| ) and SPIRE ( [Griffin &etal.|2010| ) 
perfectly matches the emission peaks of young stellar objects 
and prestellar cores, the direct progenitors of stars. It is now pos- 
sible to construct well constrained SEDs for these objects, and 
therefore determine the spatial distributions of the temperature 
and column density towards a wide range of different objects. 

Hi-GAL is the open time key project to observe the inner 
part (|/| < 60°;\b\ < 1°) of the Galactic plane with Herschel 



(Molinari & Hi-GAL Consortium 2010). The sensitivity of this 
survey will allow the detection of 10 M molecular clouds at a 
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distance of 4 kpc. A large number of diverse studies will result 
from this survey covering scientific themes such as evolved stars, 
diffuse interstellar medium, and star formation. In this Letter, 
we present the first Herschel study analyzing the temperature 
and column density structure of infrared dark clouds (IRDCs). 
These objects are reservoirs of cold, dense molecular gas seen 
in extinction in the mid-infrared against the strong background 
emission of the Galactic plane (Perault et al. 1996; Rathborne 
et al. 2006, Fig. [TJ. This gas is only slightly processed by current 
star formation activity and most likely still contains the initial 
conditions of star formation imprinted within them. Their study 
can therefore illuminate the earliest stages of star formation. In 
Sect. 2 of this paper, we describe the data, Sect. 3 focuses on the 
background properties. Section 4 presents the first results, while 
Sects. 5 and 6 summarize the discussion and conclusions of the 
paper, respectively. 



2. Data 

This analysis uses data from the Hi-GAL Science Demonstration 
Phase (SDP) which consists of two tiles of ~ 2.2° x 2.2° centered 
on / = 30° and / = 59°. For each tile, the parallel mode of obser- 
vations provides PACS 70/160 /mi and SPIRE 250/350/500 yum 
images. The SDP data w ere reduced usin g the ROMAGAL data 
reduction software (see Molinari et al. (2010) for a complete 
description). The recommended calibration factors and offsets 
were applied to the data to set the zero flux level Bern ard et al. 
( [2010] ). 
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Fig. 1. (Left): Two color image of one of the 22 IRDCs analyzed 
here. The blue is the 160 /mi emission while the red and the 
contours show the 250 /mi .(Right): Spitzer 8 /mi image of the 
same IRDC on which we have overlaid the Herschel 250 /mi 
contours. 



Fig. 2. Comparison of the original 160 /mi image (left) with 
the corresponding reconstructed background image (right). The 
intensity scale is the same for both images. The reconstructed 
background image has been smoothed to 36", the common reso- 
lution used in for the SED fitting (cf Sects. 4 and 5). The contour 
corresponds to the limit between IRDC and background derived 
from the 500 /mi data. 



In total, |Peretto & Fuller| ( |2009| ) catalogued ~ 450 Spitzer 
IRDCs in the 2 tiles of the Hi-GAL SDP fields, 75% of these 
being in the / = 30° tile. Checking each of these clouds by 
eye, we identified 80 to 90% of these IRDCs seen in emission 
with Herschel at long wavelengths, confirming that these are real 
clouds, not caused by artifacts in the Spitzer 8 /mi background. 
For this study, we selected a subsample of 22 IRDCs (see Fig. [I]), 
with sizes ranging from 40" to 180"[|]all large enough to contain 
at least one 500 /mi beam (i.e., 36") and with a column density 
peak derived from extinction such as N^ ea > 3 x 10 22 cm -2 , 
which is high enough to clearly stand out of the Hi-GAL data 
background. All images were resampled to a common pixel size 
of 2". We exclude d the W43 region bec ause of its complex back- 
ground structure ( Bally & et al.||2010 ). Our selected sample is 
clearly biased towards rather large and massive IRDCs. 



3. Background emission towards the IRDC 

The equation governing the emission I v towards an IRDC in the 
optically thin case is: 



= Z? v (7f DC )XT' RDC +/ v bg 



-back -t^° c , p 



fore 



CD 

(2) 



where B v (T d , c ) is the Planck function at the temperature of 



the IRDC averaged along the line of sight, rj, is the opac- 
ity of the IRDC at the corresponding frequency, and I v g is the 
combined background and foreground emission, which is ap- 



proximately independent of 



,-IRDC 



in the optically thin limit. 



Equation [2] shows that to calculate the opacity and temperature 
for images at a range of wavelengths, the background emission 
I v g at each wavelength needs first to be determined. 

The full images show that a significant fracti on of the tota l 
flux is from large-scale extended emission (Molinari et al. 2010). 
It is important to account for this emission to reliably map the 
IRDC emission because fluctuations in this background could 
otherwise be interpreted as structure in the IRDC. However 



1 These sizes are those derived from the extinction and given in 
Peretto & Fuller (2009). The sources are systematically larger in emis- 
sion because the Hi-GAL data probes lower column density material. 



defining the background structure is not trivial. Since the cir- 



rus noise decreases towards longer wavelengths (Gau tier et al. 
1992), we used the 500 /mi image to define the boundaries of 
each IRDC (The 350/mi data, which are very similar but have 
higherr angular resolution, were also used to confirm the bound- 
ary of each IRDC.) The pixel flux distribution at 500 /mi toward 
IRDCs shows in most cases a well defined peak at low fluxes 
and extended tail at larger fluxes corresponding to the IRDC. 
The level to separate the background and IRDC was chosen to 
be one standard deviation above the peak value, the standard de- 
viation being calculated from only those pixels below the peak. 
This procedure allows us to identify the pixels in the image as- 
sociated with the IRDC, and those pixels containing only back- 
ground emission (Fig. [2]). This procedure works well, even for 
small IRDCs, provided there are sufficient pixels to characterize 
the background emission. 

To take the background fluctuations into account, we recon- 
structed at all wavelengths background images by interpolating 
the background pixels at the position of the IRDC pixels. After 
experimenting with various interpolation methods, we decided 
to use an interpolation based on the nearest neighbors. However, 
the interpolation method has little impact on the global physi- 
cal properties since flux uncertainties are dominated by cirrus 
noise at short wavelengths (cf. next paragraph) and calibration 
uncertainties at longer wavelengths. In Fig. [2j we show the re- 
constructed background images of the IRDC at 160 /mi as well 
as the original image on the same intensity scale. The most sig- 
nificant feature of this background image is the important de- 
crease in intensity on the left part of the image (as also seen at 
250 /mi ; Fig. [T]), which is recovered well on the reconstructed 
background image. 

Uncertainties in these background images are obviously 
higher at the center of the IRDCs, where pixels are further away 
from a background pixel. To estimate this uncertainty, we com- 
puted a two point flux difference for all background pixels as a 
function of their angular separation, and this at each wavelength 
(Fig. [3]). As observed for all IRDCs, the uncertainties are much 
higher at 160 /mi than at 500 /mi , justifying our choice of choos- 
ing the latter for defining the IRDC boundaries. When the uncer- 
tainties computed this way are higher than the 20% calibration 
uncertainty, we used them on the flux measurement when per- 
forming pixel by pixel SED fitting (cf. Sect. 4). 
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Fig. 3. Evolution of the flux difference between 2 background- 
only pixels as a function of their angular separation for the IRDC 
shown in Fig. [T] (after smoothing to a common 36" resolution). 
This clearly shows that the background is far more structured at 
160 /mi than at 500 yum , justifying our choice of using the latter 
to define the boundaries of the IRDCs. 



The SED of the background images were fitted to determine 
the properties of the dust responsible for the background emis- 
sion (see Fig. [4]). With 5 data points per pixel (from 70 to 500 
yum ), we can constrain, in addition to opacity and dust tempera- 
ture, the spectral index ft of the specific dust opacity law, where 
t v oc/. Doing so, we find that the background dust tempera- 
ture varies from ~ 20 up to ~ 30 K depending on the position 
of the IRDC in the Galactic plane, in agreement with Bernard 
et al. (2010), while j3 varies from 1.6 up to 2.3, consistent with 
Boulangeretal.|(|1996|). 



4. Column density and dust temperature maps of 
IRDCs 

The reconstructed background images provide the I v g term in 
Eq.[2]for each pixel of the IRDC. The values of rf DC and rf DC 
at each pixel are then determined by minimizing the^f 2 between 
the observed fluxes and those given by Eq. |2] For this purpose, 
we used the MPFITS package (Markwardt 2009), as for the 
background SED fitting. 

The SED fitting for IRDC pixels was performed using 4 data 
points per pixel, smoothed to the same resolution of 36". The 
70 yum point was excluded since at this wavelength the optically 
thin assumption is not valid, the IRDCs are usually seen in ab- 
sorption and, therefore, the optically thin assumption is invalid 
(see Sect. 5 for a discussion about 70 yum point like sources). 
Because of the small number of data points for each fit, we 
fixed ft leaving two free parameters. Adopting the same value 
of/3 = 2 for all IRDCs, we calculated dust temperature maps 
(cf Figs. [4] and [5]). The first important feature of the map is that 
the dust temperature is non-uniform, ranging from 10 K for the 
coldest central cores, up to 22 K, consistent with the inferred 
background temperature. 

To obtain a superior angular resolution, we iterated the SED 
fitting excluding this time the 500 yum point, and using the dust 
temperature map and column density map calculated in the first 
iteration as an input parameter to the fit of the second iteration. 
Doing so, the temperature map changes little, while the column 
density map appears more structured. This is not really surpris- 
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Fig. 4. SED fit of a background-only pixel of the IRDC shown in 
Fig. [T (red) and the pixel at the column density peak (blue). The 
solid lines are the best fit to the data points. The uncertainties 
in the background are the 20% calibration ones, while for the 
IRDC we took the larger of the calibration uncertainty and the 
ones described by Fig. [3] 

ing because the temperature is predominantly constrained by the 
shortest wavelength data, toward the peak emission and on the 
exponential part of the Planck function. Given this, we directly 
used the 36" resolution temperature map in combination with 
the 250 /mi images to obtain 18" resolution column density 
maps. The resulting column density map for the IRDC shown in 
Fig. [T] is shown as contours in Fig. [5] We assumed a normaliza- 
tion of the specific dust opacity law of 0.12 cm 2 g " 1 at 250 /mi , 
adopting the thi n ice mantle dust grain models of |Qsse nkopf & 
Henning (1994) at a density n=10 5 cm -3 . Uncertainties in this 
normalization and J3 directly implies at least a factor of 2 of un- 
certainty in the column density and a couple of degrees for the 
temperatures. It is possible that J3 decreases systematically to- 
wards the i nner region of IRDCs because of dust coagulation 
(Ossenkopf & Henning 1994). However, even a decrease from 
p = 2 to ft = 1.5 increases the minimum temperature by only 
~ 2 K, which is insufficient to reverse the observed trends dis- 
cussed in the next section. 

We can already note several important points here. First, the 
strongest peaks at 250 /mi seen in Fig. [T] are not the strongest 
column density peaks. Dust temperature variations have an im- 
portant role in determining the emission. Most of the column 
density peaks are also correlated with local temperature minima, 
although there is not a one to one correlation. The 3D structure 
of the cloud is probably responsible for the apparent projection 
offsets. 



5. Discussion 

The most striking feature of our dust temperature maps is the 
non-isothermal nature of IRDCs. In Fig. [6] we plot the mass- 
weighted dust temperature versus the average H 2 column density 
where the blue symbols correspond to the 22 analyzed IRDCs, 
while the red symbols correspond to the dust temperature and H 2 
column density peaks we visually identified within the IRDCs 
(for instance 6, peaks are identified in the IRDC of Fig. [5]). 
While the median mass-weighted temperature is around 15 K, 
we see that it can be as low as 8 K in the inner part of IRDCs. 
Knowing that background temperature lies in a 20 to 30 K range, 
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Fig. 5. Dust temperature map (color scale) for the IRDC shown 
in Fig. [T] The contours represent the H 2 column density of this 
cloud, from 1 x 10 21 to 3.5 x 10 22 cm -2 . We see that the strongest 
peaks in the 250 yum image of Fig.[T]are not the strongest column 
density peaks. 



this means that the temperature decreases by 10 to 15 K from 
edge to center. Therefore, the dust temperature in dense molec- 
ular clouds is highly non-uniform. 

Temperature and its fluctuations are fundamental parameters 
for fragmentation in molecular clouds. Hydrodynamical simu- 
lations including radiative transfer (Bate 2009; Krumholz et al. 
2010 ) showed that protostellar heating could significantly reduce 
fragmentation on small scales. The picture we can draw from the 
study presented here suggests that IRDCs in the Galactic plane 
probably form from warm molecular gas with a temperature of 
~ 20 K, which cools efficiently to 10 K. On the other hand, 
in the previously mentioned simulations temperature gradients 
only develop once the first protostars form, in contrast to what 
we found here. 

The filled symbols in Fig. [6] correspond to IRDCs (blue) or 
column density peaks (red) that contain at least one 70 yum point- 
like source. We can see that two thirds of the IRDCs and one 
third of the column density peaks are currently forming stars, 
consistent with the statistics from Spitzer observations of a large 
sample of IRDCs (Peretto & Fuller, in prep). For these star- 
forming cores, the column density and temperature are probably 
overestimated and underestimated, respectively, since we did not 
take into account the data points at wavelengths below 100 /zm . 
However, since the longest wavelengths have more weight in the 
SED fitting (due to their lower uncertainties) the column density 
of the cold gas should not be affected too much. Another inter- 
esting point is that when extrapolating the linear relation found 
by |Dunham et al. (2008) between the 70 jim flux of low mass 
protostars and their bolometric luminosity to the 70 /mi sources 
we identified here, we find that their bolometric luminosity lies 
somewhere between a few 1000 L and a few 10 000 L© de- 
pending on the source and their distance (typically between 3 
and 6 kpc). This makes these sources excellent high-mass proto- 
stellar candidates. 

Finally, Fig. [6] clearly shows high column density fragments 
with no 70 /im point-like sources within them. These are inter- 
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Fig. 6. Dust temperature versus column density plot. The blue 
square symbols show the average column density and mass 
weighted temperature for the 22 IRDCs analyzed. The red cir- 
cular symbols show the temperature and column density at the 
column density peaks within the IRDCs. The filled symbols cor- 
respond to IRDCs with at least one 70yum point-like source, i.e., 
star-forming IRDC/peaks. 



esting objects for high-resolution follow-up studies since they 
might well be the birth places of massive stars, before any star 
formation activity has started. Extrapolated to the entire Spitzer 
IRDC catalog of |Peretto & Fuller | ( [2009] ), we expect to find a 
few hundred similar starless high column density cores in the 
full Hi-GAL dataset. 



6. Conclusions 

We have developed a technique to obtain reliable temperature 
and column density maps of IRDCs using the first Herschel data 
from the Hi-GAL open-time key project. We have demonstrated 
our method for a small sample of massive IRDCs by show- 
ing that they are non-isothermal but have a definite temperature 
structure, their inner regions being cooler than their edges. This 
temperature structure probably has an important impact on the 
fragmentation of the cloud in the early stages of its evolution. We 
also identified a number of high column density peaks, reaching 
values of 10 23 cm -2 , which are most likely to be the birth sites 
of massive stars. The full Hi-GAL dataset is expected to contain 
several hundreds of these sources. 
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